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Abstract 

We study SU(2) lattice gauge theory with two flavours of Dirac fermions in the 
adjoint representation. We measure the running of the coupling in the Schrodinger 
Functional (SF) scheme and find it is consistent with existing results. We discuss 
how systematic errors affect the evidence for an infrared fixed point (IRFP). We 
present the first measurement of the running of the mass in the SF scheme. The 
anomalous dimension of the chiral condensate, which is relevant for phenomenolog- 
ical applications, can be easily extracted from the running of the mass, under the 
assumption that the theory has an IRFP. At the current level of accuracy, we can 
estimate 0.05 < 7 < 0.56 at the IRFP. 



1 Introduction 



Experiments at the LHC are about to probe nature at the TeV scale, where new physics 
beyond the Standard Model (BSM) is expected to be found. The existence of a new 
strongly-interacting sector that is responsible for electroweak symmetry breaking is an 
interesting possibility. Technicolor was originally proposed thirty years ago, and strongly- 
interacting BSM has been revisited in many instances since then. Recent reviews can be 
found in Refs. [BE]. 

In order to be phenomenologically viable, technicolor theories need to obey the con- 
straints from precision measurements at LEP [21 Hj. Moreover the symmetry breaking 
needs to be communicated to the Standard Model, so that the usual low-energy physics 
is recovered. This is usually achieved in the so-called Extended Technicolor (ETC) mod- 
els by invoking some further interaction at higher energies that couples the technicolor 
sector to the Standard Model. At the TeV scale the remnants of this coupling are higher- 
dimensional operators in the effective Hamiltonian, which are suppressed by powers of the 
high energy scale, M, that characterises the extended model. Amongst these operators 
are a mass term for the Standard Model quarks, and four-fermion interactions that would 
contribute to flavour-changing neutral currents (FCNC). Thus there is a tension on the 
possible values of M: on the one hand M needs to be large so that FCNC interactions 
are suppressed, on the other hand M needs to be small enough to generate the heavier 
quark masses. In particular, the effective operator for the Standard Model quark masses 
is: _ 

C^Wf (i.i) 

where ip indicates the quark field, and $ is the field in the technicolor theory which is re- 
sponsible for electroweak symmetry breaking. In the traditional technicolor models, which 
are realised as SU(iV) gauge theories, $ = ^>^> is the chiral condensate of techniquarks. 
Let us emphasise that quark masses are defined in a given renormalisation scheme and at 
a given scale. For instance the data reported in the Particle Data Group summaries [3] 
usually refer to the quark mass in the MS scheme at 2 GeV. The coefficient that appears 
in Eq. (11.11) is the chiral condensate at the scale M: 

•M 



(f*>L=<**>| A «p / ^y(A0 

Ua f 1 



1.2) 



where 7 is the anomalous dimension of the scalar density, and A is the typical scale of 
the technicolor theory, A « 1 TeV. The chiral condensate at this scale is expected to be 
~ A 3 , and therefore the naive expectation for the quark masses is m ~ A 3 /M 2 . 
Eq. (11. 2p suggests a possible way to resolve the tension due to the large quark masses. 
If the technicolor theory is such that 7 is approximately constant (and large) over a suffi- 
ciently long range in energies, then the running above will generate a power enhancement 
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of the condensate. This scenario has been known for a long time under the name of walking 
technicolor [SI HE]- Gauge theories with a large number of fermions have been traditional 
candidates for walking theories; the fermions slow down the running of the coupling and 
can potentially lead to the required power-enhancement. More recent incarnations have 
been proposed that are constructed as SU(N) gauge theories with fermions in higher- 
dimensional representations of the colour group P QUI ED]- These theories could have a 
genuine IR fixed point (IRFP), or simply lie in its vicinity. The existence of an IRFP 
is a difficult problem to address since it requires to perform quantitative computations 
in a strongly-interacting theory. Lattice simulations can provide first-principle results 
that can help in determining the phenomenological viability of these models; numerical 
simulations of models of dynamical electroweak symmetry breaking have attracted grow- 
ing attention in recent years [121 U31 QSl US HSl HZl HH1 HSl 1201 I2H 1221 1231 1231 
EZll2Hll22lEQll^l32ll^E2EaESlEZlESlE2lSD!. A number of theories have been 
studied: SU(3) with 8, 10, 12 flavours of fermions in the fundamental representation, 
SU(3) with fermions in the sextet representation, and SU(2) with fermions in the adjoint 
representation. These studies have focused either on the spectrum of the theories, or on 
the running of the coupling computed in the Schrodinger functional (SF) scheme, finding 
some tantalising numerical evidence for IR behaviours different from what is known from 
QCD. 

Existing simulations of the Schrodinger functional have identified a possible fixed point 
in many of the above-mentioned theories by noticing a flat behaviour of the running 
coupling in this scheme over a given range of energy scales. 

In this work we consider the SU(2) theory with two flavours of adjoint fermions, and 
compute the running coupling in the SF scheme. We confirm the results obtained in 
Ref . [32J , and present a more refined analysis of the lattice data. We focus on the running 
of the mass in the SF scheme, from which we can extract the mass anomalous dimension 
that appears in Eq. (II. 2p . Current simulations are still plagued by systematic errors, 
which we examine in detail both for the coupling and the mass. These errors are the 
largest limitation to drawing strong conclusions from the lattice data. These limitations 
are common to all the studies performed so far, more extensive work is required in order 
to reach robust conclusions. Our results for the anomalous dimension of the mass provide 
crucial input for these studies that aim at exact results for non-supersymmetric gauge 
theories in the non-perturbative regime. 
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2 SF formulation 



2.1 Basic definitions 

We define the running coupling ~g 2 non-perturbatively using the Schrodinger Functional 
scheme [H], |4"2"]. This is defined on a hypercubic lattice of size L, with boundary conditions 
chosen to impose a background chromoelectric field on the system. The renormalised 
coupling is defined as a measure of the response of the system to a small change in the 
background chromoelectric field. Specifically, the spatial link matrices at t = and t = L 
are set respectively to: 

u ( x ,k)\ t=0 = exp[r)T 3 a/iL] , (2.1) 
U(x,k)\ t=L = exp[(7r-r})r 3 a/iL}, (2.2) 

with 7] = 7r/4 [13]. The fermion fields obey 

P+i) = 0, =0 at t = , (2.3) 
P-ip = 0, ^P+ = at t = L , (2.4) 

where the projectors are defined as P± = (1 ± 7o)/2. The fermion fields also satisfy 
periodic spatial boundary conditions [H]. As we mentioned above, one can readily verify 
in perturbation theory that these boundary conditions impose a constant chromoelectric 
field. 

We use the Wilson plaquette gauge action, and Wilson fermions in the adjoint repre- 
sentation, as implemented in Ref. [17] . Note that we have not improved the action, and 
therefore our results are going to be affected by 0(a) lattice artefacts. The same approach 
has been used so far for the preliminary studies of this theory in Ref. [32] . 

The coupling constant is defined as 

d^kl^-Y 1 (2.5) 



drj 

with k = — 24L 2 /a 2 sin(a 2 /L 2 (7r — 2rf)) chosen such that g 2 = g$ to leading order in per- 
turbation theory. This gives a non-perturbative definition of the coupling which depends 
on only one scale, the size of the system L. 

To measure the running of the quark mass, we calculate the pseudoscalar density 
renormalisation constant Zp. Following Ref. [33], Zp is defined by: 



Z P (L) = yWi/fp(L/2) , (2.6) 
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where f\ and fp are the correlation functions involving the boundary fermion fields ( and 





h = -l/12L 6 y d 3 ud 3 vd 3 yd 3 z(((u) l5 T a C(v)ayh 5 T a ((z)}, (2.7) 
f P (x ) = -1/12 [ d 3 yd 3 z®(x ) l5 T a il>(x )((y) l5 T a ((z)). (2.8) 



These correlators are calculated on lattices of size L, with the spatial link matrices at 
t = and t = L set to unity. 

The Schrodinger Functional boundary conditions remove the zero modes that are 
normally an obstacle to simulating at zero quark mass [JB]. This means we can run 
directly at k c . We determine k c through the PCAC mass in units of the inverse lattice 
spacing am(L/2), where 



5<So + «/aW 



and 

fA(x ) = -1/12 / d 3 yd 3 z&(x ) lo75 T^(x )ay) l5 T a ((z)}. (2.10) 



Here the lattice derivatives do and <9q are defined by dof(x) = f(x + 1) — f(x) and 
dof( x ) = f( x ) ~ f{ x ~ I); an d the correlators are calculated on lattices of size L, with 
the spatial link matrices at t — and t = L set to unity. 

We define k c by the point where am vanishes. We measure am for 5 values of k in the 
region —0.2 < am < 0.2 and use a linear interpolation in k to find an estimate of k c . The 
error on k c is estimated by the bootstrap method. 

In practice we achieve \am\ < 0.005. We check explicitly that there is no residual 
sensitivity to the small remaining quark mass by repeating some of our simulations at 
moderately small values of am ~ 0.02, for which we found no shift in ~g 2 or Zp within the 
statistical uncertainty of the measured values, so the effect of our quark mass can safely 
be neglected. 



2.2 Lattice parameters 

We have performed two sets of simulations in order to determine the running coupling 
and Zp. The parameters of the runs are summarised respectively in Tab. [U and[2J The 
values of obtained from the PCAC relation as described above. 
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L=6 


L=8 


£=12 


£=16 


P 


L=6 


L=8 


£=12 


L=16 


2.00 


0.190834 


- 


- 


- 


2.80 


0.165932 


0.165550 


0.16505 


- 


2.10 


0.186174 








3.00 


0.162320 


0.162020 


0.161636 


0.161636 


2.20 


0.182120 


0.181447 


0.1805 




3.25 


0.158505 




0.1580 




2.25 


0.180514 


0.179679 






3.50 


0.155571 


0.155361 


0.155132 


0.155132 


2.30 


0.178805 


0.178045 






3.75 


0.152803 








2.40 


0.175480 


0.174887 






4.00 


0.150822 


0.150655 






2.50 


0.172830 


0.172305 


0.17172 


0.17172 


4.50 


0.147250 


0.14720 


0.14712 


0.14712 


2.60 


0.170162 


0.169756 






8.00 


0.136500 


0.13645 


0.136415 




2.70 


0.167706 



















Table 1: Values of (3, L, k used for the determination of g 2 . The entries in the table are 
the values of k c used for each combination of f3 and L. 



Note that Zp is determined from a different set of runs at similar values of j3, L, k. 
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2.00 


0.190834 








2.80 


0.165932 


0.165550 


0.16505 




2.05 


0.188504 




0.18625 




3.00 


0.162320 


0.162020 


0.161636 


0.161636 


2.10 


0.186174 








3.25 


0.158505 




0.1580 




2.20 


0.182120 


0.181447 


0.1805 




3.50 


0.155571 


0.155361 


0.155132 


0.155132 


2.25 


0.180514 


0.179679 






3.75 


0.152803 








2.30 


0.178805 


0.178045 






4.00 


0.150822 


0.150655 


0.15051 




2.40 




0.174887 






4.50 


0.14725 


0.14720 


0.14712 


0.14712 


2.50 


0.172830 


0.172305 


0.17172 


0.17172 


8.00 


0.13650 


0.13645 


0.136415 


0.136415 


2.60 


0.170162 


0.169756 






16.0 


0.1302 


0.1302 


0.1302 


0.130375 


2.70 


0.167706 



















Table 2: Values of (3, L, k used for the determination of Zp. The entries in the table are 
the values of k c used for each combination of f3 and L. 



3 Evidence for fixed points 

Recent studies have focused on the running of the SF gauge coupling, and have highlighted 
a slow running in the lattice data for this quantity [T3J, [15], EH [32]. This is clearly different 
from the behaviour observed in QCD-like theories [131 HZ] • These results are certainly 
encouraging, but have to be interpreted with care. Lattice data can single out at best a 
range of energies over which no running is observed. However it is not possible to conclude 
from lattice data only that the plateau in the running coupling does extend to arbitrarily 
large distances, as one would expect in the presence of a genuine IRFP. On the other 
hand, if the plateau has a finite extent, i.e. if the theory seems to walk only over a finite 
range of energies, then the behaviour of the running coupling in the absence of a genuine 
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fixed point depends on the choice of the scheme, and therefore the conclusions become 
less compelling. 

Let us discuss the scheme dependence of the running coupling in more detail. The 
quantities we are interested in are the beta function and the mass anomalous dimension: 

Pj-W) = P(9), (3-1) 

l*-j-m(p) = -l(j9)™(jJ>) , ( 3 -2) 

where g, m are the running coupling and mass in a given (mass- independent) renormali- 
sation scheme. Note that 7 in Eq. (13.21) is the anomalous dimension of the scalar density, 
which appears also in Eq. (jl.2p ; 7 differs from the usual mass anomalous dimension by 
an overall sign. Both (3 and 7 can be computed in perturbation theory for small values 
of the coupling constant: 

(3(9) = -g 3 [fa + frg 2 + fat + 0(f)] , (3.3) 
j(g) = f[d + d 1 g 2 + O(t)} • (3.4) 

The coefficient /3o,Pi,do are scheme-independent; expressions for /3o,Pi for fermions in 
arbitrary representations of the gauge group have been given in Ref. [H], while for the 
first coefficient of the anomalous dimension, we have: 

d ° = (4ttV ' (3 ' 5) 
where C^-R) is the quadratic Casimir of the fermions' colour representation. In the specific 
case we are studying in this work d = 3/(47r 2 ). 

Different schemes are related by finite renormalisations; the running of the couplings in 
going from one scheme to the other is readily obtained by computing the scale dependence 
with the aid of the chain rule. Let us consider a change of scheme: 

9' = <l>(g,m/ii), (3.6) 
in' = fh!F{g : fhl '//) . (3.7) 

We impose two conditions on <p: it must be invertible, and should reduce to <f>(g) = 
g + 0(g 3 ) for small values of g. Eq. (13. 7p encodes the fact that a massless theory remains 
massless in any scheme. The picture simplifies considerably if one considers only mass- 
independent renormalisation schemes; the functions and T only depend on the coupling 
g, and one finds: 

w) = (3.8) 

i(g) = i{g)+p{g)^°%Hg)- (3-9) 
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The scheme-independence of the coefficients /3 , (3i, d can be obtained by expanding the 
functions that describe the mapping between the two schemes, <fi and T, in powers of 
~g 2 . Eqs. (13. 8p . (13.91) summarise the main features that we want to highlight here. The 
conditions we imposed on imply that ^<p{g) > 0, i.e. asymptotic freedom cannot 
be undone by a change of scheme. The existence of a fixed point is clearly scheme- 
independent: if (3{g*) = for some value g* of the coupling, then /3' has also a zero. Note 
that the value of the critical coupling changes from one scheme to the other, g'* = <fi{g*), 
however the existence of the fixed point is invariant. Similarly, the anomalous dimension 
is scheme-independent at a fixed point, since the second term in Eq. (13. 9p vanishes there. 
Moreover, if the change of scheme only involves a redefinition of the coupling, but leaves 
the mass unchanged, then the anomalous dimension does not vary. 

Unfortunately none of these conclusions holds in the absence of a fixed point. In 
particular, a flat behaviour of the running coupling over a finite range of energies can be 
obtained in any theory by a suitably-chosen change of scheme. 

It is worth stressing here another important point concerning the numerical studies 
of running couplings. There are instances where the beta function of an asymptotically 
free theory remains numerically small. This is the case of the theory considered in this 
work, namely SU(2) with 2 flavours of adjoint Dirac fermions, in the perturbative regime. 
In this case the running of the coupling is very slow from the very beginning, and this 
is independent of the possible existence of an IRFP at larger values of the coupling. As 
a consequence high numerical accuracy is needed in order to resolve a "slow" running; 
therefore numerical studies of potential IRFP need high statistics, and a robust control 
of systematics. In particular it is important to extrapolate the step-scaling functions 
computed on the lattice to the continuum limit, in order to eliminate lattice artefacts 
which could bias the analysis of the dependence of the running coupling on the scale. 
This is particularly relevant for the studies of potential IRFP, since lattice artefacts could 
more easily obscure the small running that we are trying to resolve. Some of these 
difficulties were already noted in Ref. [32J; current results, including the ones presented 
in this work, are affected by these systematics. 

More extensive simulations are therefore needed in order to remove the lattice artefacts 
by performing a controlled extrapolation of the lattice step scaling functions defined below 
in Sect. HI The scale L at which the coupling is computed and the lattice spacing a must 
be well separated. This last step is a crucial ingredient in the SF scheme, since it decouples 
the details of the lattice discretization from the running of the couplings at the scale L 
that we want to determine. Asymptotically free theories are effectively described by a 
perturbative expansion at small distances. In this regime, the degrees of freedom are the 
elementary fermions and the gauge bosons, renormalized couplings can be computed in 
perturbation theory, and different schemes can be related by perturbative calculations. 
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The evolution of the running coupling can be followed starting from this high-energy 
regime and moving towards larger distances. If the theory has an IRFP, the value of 
the running coupling approaches some finite limit g* as L is increased, i.e. the running 
coupling must lie in the interval [0, g*\. Its running can be traced from the UV regime up 
to the limiting value, which is approached from below. Larger values of g can be obtained 
in a lattice simulation; however the interpretation of these points is less transparent. One 
possibility is that the lattice theory in some region of bare parameter space lies in the basin 
of attraction of some non-trivial UV fixed point where a different continuum theory can 
be defined. The running coupling would then approach the IRFP value from above. The 
non-trivial UV fixed point is clearly difficult to identify, thereby making the extrapolation 
to the continuum limit rather tricky in this case. 

A more pragmatic approach could be to ignore the issue of the existence of a non- 
trivial UV fixed point, and simply explore the limit L/a ^> 1, assuming that the starting 
point is the lattice theory with a cutoff, and that we are only interested in the regime 
where distances are large compared to the cutoff. This interpretation is prone to sys- 
tematic errors due to potential 0(Aa) term, where A is some physical mass scale in the 
theory. These terms are not necessarily small, even if the limit a/L — > is considered. 
Moreover, the lack of a perturbative expansion prevents us defining the running coupling 
properly. The conclusion is that results for g > g* could be affected by non-universal 
lattice artefacts. 

Studies of the running couplings in the SF scheme are a useful tool to expose the 
possible existence of theories that show a conformal behaviour at large distances. However, 
the results of numerical simulations have to be interpreted with care; they are unlikely 
to provide conclusive evidence about the existence of a fixed point by themselves, but 
they can be used to check the consistency of scenarios where the long-range dynamics 
is dictated by an IRFP. A more convincing picture can emerge when these analyses are 
combined with spectral studies [13 [HI [301 HQ] , or MCRG methods [3?]. 

4 Results for the coupling 

We have measured the coupling g 2 (/3, L) for a range of /3, L. Our results are reported in 
Tab. El and plotted in Fig. [TJ it is clear that the coupling is very similar for different L/a 
at a given value of /3, and hence that it runs slowly. 

In Fig. [2] we compare our results to those obtained in Ref. [32] • Our results are directly 
comparable since we use the same action and definition of the running coupling, and it 
is reassuring to see that they agree within statistical errors. The numbers reported in 
the figure have been obtained using completely independent codes; they constitute an 
important sanity check at these early stages of simulating theories beyond QCD. 
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The running of the coupling is encoded in the step scaling function a(u, s) as 

E(u,s,a/L) = 9 2 (9o,sL/a)\_ 2{goL/a)=u , (4.1) 
a(u,s) = lim T,(u,s,a/L) , (4.2) 

a/L~>0 

as described in Ref. jl2] . The function a(u, s) is the continuum extrapolation of E(it, s,a/L) 
which is calculated at various a/L, according to the following procedure. Actual simula- 
tions have been performed at the values of f3 and L reported in Tab. [TJ 

5.0 I — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — 
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L/a=10 2/3 
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L/a=9 
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L/a=8 
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L/a=6 





_l I I l_ 



J I I l_ 



_l I I L 
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3.0 



4.0 



_L 
5.0 

P 



j i i i_ 



6.0 



7.0 



_i i i i_ 



8.0 



Figure 1: Data for the running coupling as computed from lattice simulations of the 
Schrodinger functional. Numerical simulations are performed at several values of the 
bare coupling {3, and for several lattice resolutions L/a. The points at L/a — 9, 10| are 
interpolated. 
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Table 3: Measured values of g 2 on different volumes as a function of the bare coupling (3. 
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Figure 2: The results of our numerical simulations are compared to recent results obtained 
in Ref. [32] • Different symbols correspond to different values of the lattice bare coupling 
f3, corresponding respectively to (3 — 2.5,3.0,3.5,4.5,8.0. Empty symbols correspond to 
the data obtained in this work. Full symbols correspond to the data in Ref. [32]. Symbols 
have been shifted horizontally for easier reading of the plot. 



Starting from the actual data, we interpolate quadratically in a/L to find values of 
7? 2 (/?, L) at L — 9, 10|, so that we obtain data for four steps of size s = 4/3 for L —>■ sL: 
L = 6, 8, 9, 12; sL = 8, 10|, 12, 16. Then for each L we perform an interpolation in (3 using 
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the same functional form as Ref. |31j : 



1 (3 



n 




g 2 (P,L/a) 2N 



(4.3) 



i=0 



We choose to truncate the series with the number of parameters that minimises the \ 2 
per degree of freedom. 

All the subsequent analysis is based on these interpolating functions, and does not 
make further use of the original data. Using the fitted function in Eq. (14.31) . we compute 
£(u, 4/3, a/L) at a number of points in the range u G [0.5, 3.5]. A continuum extrapolation 
is then performed in a/L using these points to give a single estimate of <j(u) = a(u,A/3). 
Example extrapolations for three values of u are shown in Fig. [3j The L = 6 data were 
found to have large 0(a) artifacts, and are not used in the continuum extrapolation. The 
L = 16 data have a large statistical error, which limits their current impact on the con- 
tinuum extrapolation. The sources of systematic uncertainty in our final results for a(u) 
are due to the interpolation in L and (3 and to the extrapolation to the continuum limit. 
Full details of the statistical and systematic error analysis are provided in Appendix lAl 

The resulting values for o(u) with statistical errors only can be seen as the black 
circles in Fig. HI The red error bars in Fig. H] also include systematic errors, but using 
only a constant continuum extrapolation. This is equivalent to the assumption that lattice 
artefacts are negligible in our data. A similar assumption has been used in Ref. [32], where 
the data at finite a/L were used directly to constrain the parameters that appear in the (3 
function of the theory. The study of the lattice step scaling function, and its continuum 
extrapolation, that we employ for this work, will ultimately allow us to obtain a full 
control over the systematic errors. 

The step scaling function encodes the same information as the (3 function. The relation 
between the two functions for a generic rescaling of lengths by a factor s is given by: 



The step scaling function can be computed at a given order in perturbation theory by 
using the analytic expression for the perturbative (3 function, and solving Eq. (14.41) for 
a(u, s). On the other hand, it can be seen directly from the definition of a(u, s) in Eq. (14.21) 
that an IRFP corresponds to a(u, s) = u. 

Our current values for the step scaling function are consistent with a fixed point in the 
region ~g 2 ~ 2.0 — 3.2, as reported in Ref. [32]. Further simulation at higher ~g 2 is limited 
by the bulk transition observed in Ref. [HI [30] at (3 ~ 2.0. 

The errors from also including the linear continuum extrapolation are much larger and 
mask any evidence for a fixed point, as shown in Fig. This should be a conservative 
estimate of the total uncertainty on a(u), which is dominated by systematic errors. 
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Figure 3: Results for the lattice step-scaling function £(4/3, u, a/L). The dashed lines 
represent the initial value of u. The point at x = yields the value of cr(u), i.e. the 
extrapolation of £ to the continuum limit. The error bar shows the difference between 
constant and linear extrapolation functions, and gives an estimate of the systematic error 
in the extrapolation as discussed in the text. 
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Figure 4: The relative step-scaling function u{u) / u obtained after extrapolating the lattice 
data to the continuum limit. The black circles have a statistical error only. The red error 
bars also include systematic errors, but using only a constant continuum extrapolation 
(i.e. ignoring lattice artifacts). Note that a fixed point is identified by the condition 
a{u)/u = 1. 
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Figure 5: The relative step-scaling function a{u)/u obtained after extrapolating the lattice 
data to the continuum limit. The black circles have a statistical error only, the red error 
bars include systematic errors but using only a constant continuum extrapolation, and 
the grey error bars give an idea of the total error by including both constant and linear 
continuum extrapolations. 



5 Running mass 

The running of the fermion mass is determined by the scale-dependence of the renormal- 
isation constant for the pseudoscalar fermion bilinear Zp defined in Eq. (I2.6p . Note that 
Zp is both scheme and scale dependent. The same step scaling technique described for 
the gauge coupling can be used to follow the nonperturbative evolution of the fermion 
mass in the SF scheme. In this work, we follow closely the procedure outlined in Ref. [15] . 

We have measured the pseudoscalar density renormalisation constant Z P (j3, L) for a 
range of f3, L. Our results are reported in Tab. HI and plotted in Fig. [6], where we see that 
there is a clear trend in Zp as a function of L at all values of f3. 

The lattice step scaling function for the mass is defined as: 

Z P (g ,sL/a) 



T,p(u, s, a/L) 



Zp(g , L/a) 



; (5.i) 

g 2 (L)=u 
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the mass step scaling function in the continuum limit, ap(u, s), is given by: 



<ip(u, s) = lim £p(w, s, a/L) . 

a— >0 



(5.2) 
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Figure 6: Data for the renormalisation constant Zp as computed from lattice simulations 
of the Schrodinger functional. Numerical simulations are performed at several values of 
the bare coupling (3, and for several lattice resolutions L/a. The points at L/a = 9, 10| 
are interpolated. 
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16.0 


0.9146(4) 
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0.9005(5) 


0.8887(15) 


2.70 


0.5228(6) 



















Table 4: Measured values of Zp on different volumes as a function of the bare coupling 
P. 
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The method for calculating o>(u) = <Tp(u, 4/3) is similar to that outlined in Sec. H]for 
calculating a(u). Interpolation in f3 is accomplished using a function of the form: 



M/W«) = E C *(£) 



(5.3) 



Full details of the procedure are given in Appendix [Bl Again the errors are dominated 
by systematics, in particular the choice of continuum extrapolation function. In Fig. [7] 
we see that, unlike g 2 , Zp has a significant variation with a/L that is fit well by a linear 
continuum extrapolation. The constant extrapolation is only used to quantify the errors 
in extrapolation. 
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Figure 7: Results for the lattice step-scaling function £p(4/3, u, a/L). The point at x = 
yields the value of ap(u), i.e. the extrapolation of Sp to the continuum limit. The error 
bar shows the difference between constant and linear extrapolation functions, and gives 
an estimate of the systematic error in the extrapolation as discussed in the text. 



16 



1.04 



1.02 
1.00 
0.98 
_ 0.96 

S 0.94 

O 

0.92 
0.90 
0.88 
0.86 
0.84 



1-loop 

Statistical Error 



j i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i_ 







0.5 



1.5 



2.5 



3.5 



Figure 8: The step-scaling function for the running mass ap(u), using a linear continuum 
extrapolation. The black circles have a statistical error only, the red error bars include 
systematic errors using a linear continuum extrapolation. The grey error bars come from 
also including a constant extrapolation of the two points closest to the continuum, and 
give an idea of the systematic error in the continuum extrapolation. 



Using the fact that <7p(u, s) = fn(fi) /m(/i/ s) for fi = 1/L, we can perform an iterative 
step scaling of the coupling and the mass to determine the running of the mass with scale. 
However, since we observe no running of the coupling within errors this is not particularly 
interesting. 

The mass step scaling function is related to the mass anomalous dimension (see e.g. 
Ref. [28]): 



Op (it 



a. 



a(u) 



W)/(2/3 )) 



exp 



dx 



l{x) d 



(3{x) (3 x 



(5.4) 



We find good agreement with the 1-loop perturbative prediction, as shown in Fig. [HJ 

In the vicinity of an IRFP the relation between op and 7 simplifies. Denoting by 7* 
the value of the anomalous dimension at the IRFP, we obtain: 



ra W«) dm 

m(» m 



» /s dq 
Q 



7 — 



(5.5) 
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and hence: 

log \a P (s,u)\ = -7* logs. (5.6) 
We can therefore define an estimator 

\og\a P (u,s)\ 
log \s\ 

which yields the value of the anomalous dimension at the fixed point. Away from the 
fixed point 7 will deviate from the anomalous dimension, with the discrepancy becoming 
larger as the anomalous dimension develops a sizeable dependence on the energy scale. 

We plot the estimator 7 in Fig. Again the error bars come from evaluating the above 
expression using the extremal values of <Jp{u) at each u. We see that the actual value 
of 7 is rather small over the range of interest. In particular at Tj 2 = 2.2, the benchmark 
value for the IRFP tentatively found in Ref. [32J, we have 7 = 0.116^28 using just the 
linear continuum extrapolation, and 7 = 0.116^28 if we include the constant continuum 
extrapolation as well. In the presence of an IRFP 7 yields the value of the anomalous 
dimension, and therefore the values above can be used to bound the possible values of 7*. 
The results of Ref. [32] suggest the IRFP is in the range Tj 2 = 2.0 — 3.2; at the extremes of 
this range we find 7* = 0.086^ and 0.41+^ using just the linear continuum extrapolation, 
and 7* = 0.086^Jq 5 and O.41I33 including the constant continuum extrapolation. Over 
the entire range of couplings consistent with an IRFP, 7* is constrained to lie in the 
range 0.05 < 7* < 0.56, even with our more conservative assessment of the continuum 
extrapolation errors. 
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Figure 9: The mass anomalous dimension estimator 7(w). The dashed line shows the 
1-loop perturbative result, the black circles have a statistical error only, and the red error 
bars include systematic errors using a linear continuum extrapolation. The grey error 
bars also include a constant extrapolation of the two points closest to the continuum, 
giving an idea of the systematic error involved in the continuum extrapolation. 



6 Conclusions 

In this paper we have presented results for the running of the Schrodinger Functional 
coupling ~g 2 and the mass anomalous dimension 7. 

Turning first to the running of the coupling, our results are completely consistent with 
those of Ref. [32] • Our statistical errors are larger; however, we have carried out our 
analysis in a way that aims at disentangling clearly the scale dependence from the lattice 
artefacts. Our analysis can be systematically improved as more extensive studies are 
performed, and will ultimately allow us to take the continuum limit with full control over 
the resulting systematic errors. Our results appear to show a slowing in the running of the 
coupling above g 2 = 2 or so, and are consistent with the presence of a fixed point where 
the running stops at somewhat higher ~g 2 . This is consistent with the analysis of Ref. [32] . 
However, once we include the systematic errors from the continuum extrapolation we find 
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that our results no longer give any evidence for a fixed point. The fundamental reason for 
this is that the running of the coupling is very slow in this theory and so great accuracy 
is needed, in particular near a possible fixed point. 

By contrast, we find that the behaviour of the anomalous dimension 7 is much easier to 
establish. The systematic errors from the continuum extrapolation are much smaller than 
the signal, and we find a moderate anomalous dimension, close to the 1-loop perturbative 
prediction, throughout the range of j3 explored. In particular, in the range g 2 = 2.0 — 3.2, 
where there may be an infrared fixed point, we find 0.05 < 7 < 0.56. These values are 
much smaller than those required for phenomenology, which are typically of order 1-2. 
Such large values of 7 are clearly inconsistent with our results. The anomalous dimension 
at the fixed point can be computed analytically using the all-order beta function proposed 
in Ref . . The result can be expressed as a function of group-theoretical factors only. 
Using the conventions described in the Appendix of Ref. [Hj for these group-theoretical 
factors, the result in Ref. [49] yields 7 = 3/4, which is not too far from the bound we 
quote above. Given the uncertainty in the exact value of the SF coupling at the fixed 
point, 7 = 3/4 is not strongly excluded by our current data. More precise investigations 
are needed to clarify this point. 

The anomalous dimension is easier to measure than the beta function in candidate 
walking technicolor theories since it is expected to be different from zero, so we are 
measuring the difference of two quantities that are significantly different, say Zp at L = 8 
and L = 12. By contrast for the running of the coupling we must measure the difference 
of two quantities that are almost the same, since the beta function is expected to be small. 
Furthermore, the anomalous dimension is crucial for phenomenology; if it is not large then 
the presence or absence of walking behaviour becomes academic. Hence the implications 
of our measurement of 7 for the phenomenology of minimal walking technicolor call for a 
more precise study. 

Our conclusion that 7 is not large is unlikely to be affected by using larger lattices. 
One can see this by considering the continuum extrapolations in Fig. [71 For 7 to reach, 
say, 1 in the continuum limit, we would need Sp to be 3/4 = 0.75 at a/L = 0. However 
we see that the dependence on a/L is much too small for this to be possible, and indeed 
is in the wrong direction. Only a very unlikely conspiracy of lattice artifacts would make 
it possible for Ep to be as small as 0.75 in the continuum limit. On the other hand the 
value of g corresponding to the IRFP is currently not known with sufficient accuracy. 

The results presented here are the first computation of the anomalous dimension at 
a putative fixed point; the systematic errors need to be reduced to make our conclusions 
more robust. In particular, using larger lattices would give results at smaller a/L and 
hence make the continuum extrapolations more accurate. It may also be necessary to 
use an improved action in the long term to achieve the precision required to show the 
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existence of an IRFP or of walking behaviour. However, as described above, this is very 
unlikely to affect our phenomenologically most important result, namely that 7 is not 
large. Recent results in Ref. [50] suggest that the anomalous dimension can be computed 
using finite-size scaling techniques. A comparison of different techniques will improve the 
determination of the anomalous dimension. 
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A Coupling error analysis 

We directly measure the Schrodinger Functional coupling ~g 2 and perform multiple stages 
of interpolation and extrapolation to extract the continuum step scaling function cr(u) = 
<7(u,4/3). 





L/a 
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lOf 


12 


CO 


1.113 ± 0.057 


0.967 ± 0.050 


1.010 ±0.001 


0.987 ± 0.003 


0.988 ± 0.024 


Cl 


-0.560 ±0.206 


-0.064 ±0.215 


-0.259 ± 0.001 


-0.216 ± 0.006 


-0.226 ± 0.055 


C2 


0.130 ± 0.216 


-0.307 ± 0.328 




-0.022 ± 0.003 


-0.016 ± 0.028 


C3 


0.366 ± 0.125 


0.221 ± 0.211 








C4 


-0.136 ±0.196 


-0.059 ± 0.048 








C5 


-0.364 ±0.234 










C6 


0.298 ± 0.127 










C7 


-0.064 ± 0.024 










x 1 

dof 


2.85 


2.42 


1.73 


3.45 


3.37 


do} 


8 


7 


4 


3 


4 



Table 5: Interpolation best fit parameters for g 2 . 
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In order to estimate our errors for each of these stages we perform multiple bootstraps 
of the data. The full procedure to get a single estimate of cr(u) can be summarised as 
follows: 



• Generate Nb x N a bootstrapped ensembles of the data and extract mean and error 
for each. 

• For each bootstrap, interpolate in a/L to find values at L = 9, 10|. 

• From each set of N a of these find the mean and standard deviation, to give Nj, 
interpolated data points with error bars. 

• For each of the Nj, bootstraps do a non-linear least squares fit for g 2 (/3, L) interpo- 
lation functions in (3, an example is shown in Fig. [10J 

• Use these functions to find estimates of H(a/L,u) for L = 8,9, and from this 
extract a mean and error for each a/L. 

• Perform a single weighted continuum extrapolation in a/L using these points to give 
a(u). 

This process is repeated N m times, bringing the total number of bootstrap replicas of 
the data to N a x ]V& x N m . This gives N m estimates of cr(u), from which a mean and 
1-sigma confidence interval is extracted. 
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Figure 10: Example of an interpolation function for L 
compared with measured g 2 data points. 



I, with a ±<r confidence interval, 



However, the systematic errors that result from varying the number of parameters 
in the interpolation functions or the continuum extrapolation functions are significantly 
larger than the statistical errors for the optimal set of parameters. 

In order to quantify this, we repeated the entire bootstrapped process of calculating 
<y(u) with a range of different interpolation and extrapolation functions, each of which 
gives an estimates for a(u), with a statistical error. 

Specifically, we included two different choices for the number of parameters in the 
interpolating functions at each L. We kept the best fit, outlined in Tab. Eland added the 
function with the second lowest \ 2 P er degree of freedom as shown in Tab. The error 
in the continuum extrapolation was estimated by including both constant and linear 
extrapolation functions. All possible combinations of these functions gave us a set of 
2 5 = 32 values for a(u), each with a statistical error, which spanned the range of the 
systematic variation. 

For each value of u the resulting extremal values of cr(u) were used as upper and lower 
bounds on the central value. 
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9 2 
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1.010 ±0.001 
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Cl 


-0.560 ±0.206 


-0.064 ±0.215 


-0.259 ± 0.001 
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-0.226 ± 0.055 
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-0.016 ± 0.028 
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0.366 ± 0.125 


0.221 ± 0.211 
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-0.136 ±0.196 


-0.059 ± 0.048 
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-0.364 ±0.234 
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0.298 ± 0.127 
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do} 
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3.45 


3.37 


do} 
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4 



Table 6: Interpolation next-best fit parameters for g 2 . 

B Mass error analysis 

The mass error analysis follows the same procedure as outlined in Appendix [A] with ~g 2 
replaced by Zp. The function used to interpolate Zp in (3 is given in Eq. 15. 3[ and an 
example fit is shown in Fig. [TTJ The q giving the smallest reduced y 2 are given in Tab. [7] 
and those with the second smallest in Tab. [HJ 

In addition, Zp converges faster than ~g 2 and we have better 16 4 data so we can use 
3 points in our continuum extrapolations. Again the L = 6 data were found to have 
large 0(a) artifacts so are not used in the continuum extrapolation, and for the constant 
extrapolation only the two points closest to the continuum limit are used. The fits for 
both ~g 2 and Zp are required to determine ap(u), so independently varying the choice of 
the number of parameters for these now gives 2 10 = 1024 values for <jp(u), each with a 
statistical error. 
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0.58 ± 0.30 


0.93 ± 0.09 


1.02 ±0.01 


1.00 ± 0.01 
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7.64 ±6.85 


-0.43 ± 1.74 


-2.17 ± 0.10 


-1.76 ±0.01 


-1.98 ±0.08 


-1.99 ±0.09 


C2 


-78.87 ± 60.50 
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4.70 ±0.54 
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Table 7: Interpolation best fit parameters for Zp. 
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-105.52 ± 15.92 
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Table 8: Interpolation next-best fit parameters for Zp. 
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Figure 11: Example of an interpolation function for L = 8, with a ±cr confidence interval, 
compared with measured Z P data points. 
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